AN EXAMINATION OF THE EFFECT OF ATTACK 
VELOCITY ON THE OUTCOME OF LANCHESTER-TYPE 
ENGAGEMENTS WITH RANGE DEPENDENT KILL-RATES 



By 



James Francis Lloyd 




A ppxovcd public xclcasc; diAtxibiution untimite.d. 



r ) 3 ? 3 ‘j j 



An Examination of the Effect of Attack 



Velocity on the Outcome of Lanchester-Type 
Engagements with Range Dependent Kill-Rate 



by 



James Francis^Lloyd , Jr. 
Captain, United Stages Marine Corps 
B.S., United States Naval Academy, 1963 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN OPERATIONS RESEARCH 



from the 

NAVAL POSTGRADUATE SCHOOL 
March 1971 



library 

UAVAL POSTGRADUATE SCHOOL 
M02JTEREY , CALIF. 93940 



ABSTRACT 

This thesis examines the effect of attack velocity on 
the outcome of Lanchester- type engagements between forces 
with range dependent kill-rates. Range dependent (linear 
and quadratic) kill-rates are considered, and analytic 
solutions to Lanchester- type equations are utilized in 
this study. 

By varying the attack velocity, the effects on ter- 
minal force strengths are investigated for the case when 
an attacking force has the initial fighting strength su- 
periority, and for the case when a defending force has the 
initial fighting strength superiority. 
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I. INTRODUCTION 



S. Bonder and colleagues [8] have described the neces- 
sity for military planners to evaluate, both comparatively 
and individually, the effectiveness and costs of proposed 
weapon systems and force structures. They have pointed out 
that Monte Carlo simulations or wargames have been used ex- 
tensively in situations where adequate descriptive theory 
has not existed. The shortcomings of the Monte Carlo simu- 
lations and wargames in terms of cost and degree of diffi- 
culty of analysis are noted, as are the advantages and 
disadvantages of analytic models. One of the disadvantages 
of analytic models is the limited number of them available. 
The best known of the analytic models which are useful in 
weapon system and force structure evaluation are the 
Lanchester theories of combat. 

Lanchester equations are sets of simultaneous differ- 
ential equations which describe deterministically the at- 
trition of two opposing homogeneous forces [2]. The more 
well known of the Lanchester equations are those which 
describe combat between two forces employing aimed fire 
against each other; they are: 
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where : 



a = constant rate at which a single y unit destroys a 
single x unit. 

3 = constant rate at which a single x unit destroys a 
single y unit. 

x(t) = number of X force survivors at time t. 
y(t) = number of Y force survivors at time t. 

The quantities a and 3 are known as the attrition- rate co- 
efficients, or kill-rates. The state solution of (1) and 
(2), i.e., the solution with the independent variable, time, 
removed, is: 

3(X 2 - x 2 ) = a (Y 2 - y 2 ) , (3) 

where X and Y are the initial force strengths. This is re- 
ferred to as the Lanchester square law. If the forces are 
evenly matched during an engagement, x and y, the number of 
survivors, approach zero together. Equation (3) then be- 
comes : 

3X 2 = otY 2 , (4) 

which is the parity condition between forces. It is the 
condition for a draw. 

The above Lanchester- type equations, (1) and (2), are 
based upon the following assumptions: 

1. Opposing forces consist of homogeneous units. 

2. All units of each force are within range of all 
weapons of all units of the opposing force. 

3. Each force has perfect intelligence; i.e., the 
exact location of each unit of the opposing force 
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is known; and when an opposing unit is destroyed, 
fire is immediately shifted to a surviving unit. 

4. Fire is distributed uniformly over all surviving 
targets . 

In light of these necessary assumptions, several deficiencies 
in the model exist, prohibiting the use of Lanchester equa- 
tions as effective planning tools. Only problems concern- 
ing forces composed of homogeneous units can be solved 
generally. It is sometimes difficult to theoretically pre- 
dict the attrition- rate coefficients of proposed weapon 
systems. The effects of mobility on the weapon system and 
on the attrition- rate coefficient have not been generally 
incorporated into the Lanchester model [8], 

Weiss [9] presented an extension of Lanchester theory 
of combat which included the relative movement of forces, 
and thereby allowed time and space to be "traded" for cas- 
ualties. This extension was an attempt to incorporate the 
fact that the attrition- rate coefficients are dependent on 
force separation. One of Weiss's assumptions, however, was 
that during an engagement, force separation reached an 
equilibrium point, and that a force would advance or re- 
treat based on a comparison of the actual casualty rate with 
some predetermined tolerable casualty rate. Thus, if a 
force were receiving casualties at a rate higher than the 
predetermined acceptable rate, that force would retreat; 
likewise, the other force would advance to maintain its 
equilibrium casualty rate. Weiss further assumed that the 
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time required for the forces to respond to fluctuations in 
their casualty rates was small compared to the time re- 
quired for the forces to close. Therefore, there was no 
change in force separation, and hence, no change in the 
attrition-rate coefficients. Movement of forces was in- 
cluded in the model, but the effect of force separation on 
the attrition- rate coefficients was disregarded. 

Bonder [1], [3], [8], has further developed Weiss's 

extensions of Lanchester equations to investigate the ef- 
fects of mobility and range dependencies of weapon systems. 

He has done this by formulating a model which considers 
mobility and the influence of range on the attrition-rate 
coefficients. The attrition equations of this model are 
given by: 

gf = -a(r)y(t) (5) 

= _ B(r)x(t) , (6) 

where the terms are defined as for equations (1) and (2), 
with the exception that the attrition- rate coefficients are 
no longer constant, but now vary with the range between 
forces . 

Three forms of the attrition- rate coefficient were 
considered by Bonder. They were: 

a(r)=k a (R e - r ) , 8 (r) =k g (R e -r) (linear form) (7) 

ct (r) =a 0 (1-r/R ) 2 , B (r) = 8 0 (1 - r/R ) 2 (quadratic form) (8) 

a(r)=a 0 /2[l+cos('rTr/R ) ] , 3 (r) = Bo / 2 [1 + cos (nr/R ) ] (cosine 

e e form) 

(9) 
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where : 

r = force separation, 

R = effective range of both forces, 

k a , = constants, 

ao , 3o = attrition-rate coefficients at r = 0. 

Mobility was incorporated into the model by transforming 
equations (5) and (6) into second order differential equa- 
tions with r as the independent variable [1]. By assuming 

that the effective range, R , was the same for both forces, 

c 

and the attrition- rate coefficients were of the same form, 
Bonder was able to obtain an analytic solution for the 
average force strengths, x(r) and y(r), for decreasing values 
of force separation [1] , [3] , [8] . 

Although Bonder's model is much more appropriate than 
the classical Lanchester theory for the analysis of mobile 
weapon systems, its use is delimited by the assumptions 
necessary to obtain the solution for x(r) and y(r) . Taylor 
[6], using Bonder's model, employed a different mathematical 
approach in obtaining a solution to (5) and (6) . The solu- 
tion he obtained is much more general than that obtained by 
Bonder; the assumptions necessary for Bonder's solution are 
not required for Taylor's solution. Taylor's solution is in 
the form of an infinite series, with time as the independent 
variable. Taylor has previously shown [7] that it is incon- 
sequential if either force separation or time is considered 
as the independent variable, one form being derived from 
the other. Thus, a general solution to Lanchester equations 
is now available for those situations in which the 
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opposing forces have different effective ranges and/or dif- 
ferent range dependencies of attrition-rate coefficients. 

Of particular interest is Taylor's solution to the 
Lanchester equations when the forces have attrition- rate 
coefficients of the following form: 



a(t)=k t 
v a 



m 



n 



( 10 ) 

( 11 ) 



b(t)=k b t 

The solution to the attrition equations 

jf = -a(t)y(t) 

= -b ( t) X ( t) 

with a(t) and b(t) given by (10) and (11) respectively is 
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where 



(m+n+2) 

2 



m+1 

(m+n+2) ’ 



p+q = 1 , 



s = 



P = 



(14) 



and I^(x) is the modified Bessel function of the first kind 
and order v [6] . The solution is now in terms of tabulated 
functions . 

Bonder's model and Taylor's solution for average force 
strengths given by (12) and (13) were used to extend Bonder's 
study [3], [8], of the effect of mobility on the outcome of 
a Lanchester type engagement when both forces had the same 
effective range. The attrition- rate coefficients were de- 
pendent on the force separation, but not necessarily in the 
same manner. 
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II. METHOD OF STUDY 



In investigating the effects of mobility on the out- 
come of a Lanchester- type engagement, the following ques- 
tions were considered: 

1. How were the terminal force strengths affected by 
the attack velocity? 

2. How did the form of the attrition- rate coefficient 
contribute to the effects of attack velocity? ' 

A. SCENARIO 

The following scenario was used in an attempt to answer 
the above questions. Two homogeneous forces, an X force 
and a Y force, were engaged in combat. The X force was the 
attacking force, while the Y force was defending a fixed 
position. Both forces had attrition- rate coefficients which, 
in most cases, varied with force separation. This force 
separation closed at a constant velocity. The range at 
which the battle commenced was the effective range for both 
forces. The attrition- rate coefficients used were the aver- 
age values for each value of force separation, i.e., sto- 
chastic variations in attrition rates at a specific range 
were disregarded [8] . The attrition rates used were of the 
following form: 



a(r) = a„(l-r/R e ) m 



(15) 



3 (r) = B 0 (l-r/R e ) n . 



(16) 
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Equations (15) and (16) can be shown to be equivalent to (10) 
and (11) be letting 
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where R q i s the range at which the battle begins (R =R ) , 
and V is the attack velocity. The force attrition for this 
scenario is given by (5) and (6) . The necessary assumptions 
inherent in the Lanchester theory of combat given above are 
applicable in this situation. 

B. DESCRIPTION OF MODEL 

The model used was based on Taylor's solution, given in 
(12) and (13) . Using the relationships that 
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a solution with range as the independent variable was ob- 
tained. This solution is as follows: 
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where p, s, and q are as previously defined in (14). Again, 
this solution is in terms of tabulated functions. There are, 
however, few tables of modified Bessel functions of the 
first kind of fractional order. These tables are available 
in reference [5] . The modified Bessel function of the first 
kind of order p has the power series expansion [4] 
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Using this form of the Bessel function, equations (17) and 
(18) can be rewritten as 
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A computer program which evaluated equations (19) and (20) 
was used to obtain numerical results. 



C. INPUT PARAMETERS 

Two sets of numerical results were computed using this 
model. In one case, Case I, the parameters were selected 
such that the X force had a larger initial strength than the 
Y force, but the Y force had greater fighting strength in 
the classical Lanchester sense, (i.e., apY 2 > gpX 2 ) . In the 
other case, Case II, the X force had both an initial strength 
advantage and a fighting strength advantage (gpX 2 > apY 2 ). 

The parameter values used in each of these cases are shown 
in Table I. In obtaining numerical results in each case, 
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some of the parameters were held constant, while other pa- 
rameters were varied. The parameters which were held con- 
stant were the initial force strengths, the values of ao 
and $o , and the effective range; the parameters which were 
varied were the attack velocity, and the exponents in the 
attrition- rate coefficients, m and n. For each combination 
of m and n, (19) and (20) were evaluated for values of at- 
tack velocity varying from 0.5 meters per second to 25 meters 
per second. 

The values of ao and Bo were selected arbitrarily; none 
of the results of studies of the attrition-rate coefficient 
[1], [8], were applied in obtaining these values. The values 

selected are, however, the same order of magnitude as those 
used by Bonder in his work. Figure 1 shows how the attrition - 
rate coefficients vary with range for different values of m 
and n . 

D. USE OF THE MODEL 

1 . Computer Program 

Numerical results for the model were obtained from 
a computer program written in FORTRAN IV language; this pro- 
gram is shown on pages 40 and 41 . Equations (19) and (20) 
are particularly amenable to computer based evaluation. The 
infinite series terms in both equations were evaluated using 
recursive relationships. In describing the programming 
technique used in computing these infinite series terms, 
only equation (19) will be discussed. An analagous technique 
was used for the evaluation of the terms in (20) . 
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Equation (19) can be rewritten as follows: 
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for k > 1 . 

Thus the terms for Si and S 2 were obtained from the recur- 
sive relationships for and T 2 J,. In the program, Si and 

S 2 were computed as follows: 
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where N is the smallest integer such that 1\^ < 1x10 

The average force strength for both forces was com- 
puted for decreasing incremental values of force separation. 
The program terminated when either the X force or the Y 
force was completely annihilated, or when the X force over- 
ran the Y force position, (i.e., r=0) . 

2 . Computer Program Output 

The computer program gave the average X and Y force 
strengths for every 20 meters of decreasing force separa- 
tion. From this information, the values of force strengths 
and force separations at the end of the engagement were ob- 
tained. The values of terminal force strength and final 
positions of forces were obtained for each combination of 
values of m and n, and for each value of attack velocity 



considered . 



III. DISCUSSION OF OUTPUT AND RESULTS 



Figures 2 through 10 show the effects of attack veloc- 
ity on the terminal force strengths for each combination of 
parameters considered. The graphs readily indicate that the 
effect of attack speed on the outcome of the engagement is 
dependent on the parameter values used, particularly the 
values of m and n, which determine the form of the attri- 
tion-rate coefficients. 

The classical Lanchester results are shown for Cases 1 
and II in Figure 2. When m=n=0, the attrition-rate coeffi- 
cients are constant throughout the engagement. Note that 
in Case I, i.e., when the Y force has greater fighting 
strength, the X force was annihilated before it reached Y's 
position when it attacked at speeds of less than 20 meters 
per second. At higher attack speeds, the X force was able 
to reach Y's position, but Y had numerical, and hence, 
fighting superiority. In Case II, with m=n=0, the X force 
dominated the engagement throughout. At attack speeds of 
less than 16 meters per second, the Y force was completely 
destroyed before the X force reached its position. At 
speeds greater than 16 meters per second, there were Y 
_ force survivors when the X force reached the defended posi- 
tion. As the attack speed increased, the number of X and 
Y force survivors increased, but the X force maintained 
numerical and fighting superiority. Thus, an increase in 
attack speed caused an increase in terminal strength in both 
Case I and Case II. The maximum velocity for which one side 
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was annihilated corresponds to the time required for anni- 
hilation, as obtained from the time solution of equations 
(1) and (?.) [2]. Thus, increasing the attack velocity 

merely shortens the length of the engagement. Since the 
Lanchester square law holds in this case, a decrease in the 
duration of the battle results in an increase in the number 
of survivors of both forces. Thus, by attacking at a very 
high attack speed, a force with initial numerical superior- 
ity can obtain terminal numerical superiority. 

When m=n=l, both attrition- rate coefficients vary lin- 
early with force separation. Figure 3 illustrates how the 
outcome of the engagement in Case I and Case II was affected 
by the attack velocity. In Case I, the X force was annihi- 
lated before it was able to reach Y's position when it at- 
tacked at speeds less than 10 meters per second. At veloc- 
ities greater than 10 meters per second, the number of X 
force survivors increased much faster than the number of 
Y force survivors. Finally, at velocities greater than ap- 
proximately 24 meters per second, the number of X survivors 
exceeded the number of Y survivors, resulting in X force 
numerical superiority. In Case II, an increase in attack 
velocity resulted in an increase in the number of survivors 
for both forces, and a prolongation of the engagement. The 
X force maintained its terminal numerical and fighting su- 
periority for all values of attack velocity. The situation 
when m=n=2, as shown in Figure 4, is similar to the case 
when m=n=l . In Case I, the X force achieved terminal numerical 
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superiority by attacking at a velocity greater than approxi- 
mately 16 meters per second. In Case II, an increase in 
attack velocity caused an increase in the number of survivors 
for both forces, but the X force maintained its terminal 
strength advantage at all values of attack speed. 

Thus, when both attrition- rate coefficients vary in a 
linear or quadratic manner, an attacking force can compensate 
for its lack of fighting strength superiority by attacking at 
a higher velocity. These results have been obtained by 
Bonder [3] . When the attacking force has the advantage in 
terms of fighting strength, it can increase its number of 
survivors by increasing its attack speed. This has the pos- 
sibly adverse effect of causing an increase in the number 
of enemy survivors at the objective; the attacking force, 
however, Is still the superior force, regardless of the at- 
tack speed. It would thus appear that the attacking force 
commander would have a choice of how to attack the defended 
position. A relatively slow attack would destroy the enemy 
prior to the attacker's arrival at the position, but would 
result in a larger number of casualties in the attacking 
force. Alternatively, a faster rate of attack would result 
in fewer friendly casualties, and more enemy survivors at 
the objective. 

The increase in the number of suriviors of both forces 
is a direct consequence of the Lanchester square law; the 
square law holds when both linear and quadratic kill rates 
are considered. A result of the square law is that the two 
forces will always have the same strength values at a specific 
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time during a battle, regardless of the duration of the 
battle. Since an increase in attack velocity is equivalent 
to a decrease in the length of the battle , the number of 
survivors of both forces will increase. Thus, it is pos- 
sible for a force with initial numerical superiority, but 
with firepower inferiority, to reach a defended position with 
numerical superiority by attacking at a higher speed. The 
reduction in the duration of the engagement prevents the 
defending force, with its superior firepower, from reducing 
the attacking force strength to an inferior level. A reduc- 
tion in the duration of the battle is equivalent to a trans- 
formation of the time scale of the attrition process. 

Figures 5 through 10 illustrate the outcome of the 
engagement when the opposing forces do not have attrition- 
rate coefficients of the same form. Figures 5 through 8 
represent the cases when one of the forces has a constant 
attrition-rate coefficient (m or n = 0) , and the other force 
has a linear or quadratic coefficient. As is readily obvious, 
when m=0, the attack velocity had no influence on the outcome 
of the engagement; the X force was annihilated in every in- 
stance except in Case II with m=0 and n=l. Then, as can be 
seen in Figure 5, there was a slight increase in the X force 
terminal strength at a velocity of about 25 meters per second. 
This increase was insignificant when compared with the Y force 
terminal strength. 

When n=0, i.e., the X force attrition-rate coefficient 
is constant, the attack speed did influence the outcome of 
the engagement. Figures 7 and 8 depict these cases. The X 
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force had terminal numerical superiority in each of the 
engagements, and in each case, slower attack speeds resulted 
in more X force survivors. This is reasonable, since, as 
can be seen in Figure 1, the X force attrition-rate coeffi- 
cient dominates the Y force coefficient at greater ranges 
between forces. Thus, at slower attack speeds, more time 
is spent at these greater ranges where the X force has the 
larger kill-rate. From this result, it is seen that it is 
advantageous to keep the enemy force in the best "kill zone." 

A more practical, and common situation is when one force 
has a linear attrition- rate coefficient and the other force 
has a quadratic coefficient. Figures 9 and 10 illustrate 
these situations. In both Case I and Case II, with m=l and 
n=2, the Y force had a large terminal numerical advantage 
at the slower attack speeds. However, the X force was able 
to reduce this superiority by attacking at a higher velocity. 
In Case I , the X force reduced its casualties by increasing 
its speed, even though it was never able to achieve a num- 
erical advantage. However, in Case II, the increased attack 
velocity did result in X force numerical superiority. In 
both of these cases, the X force attrition- rate coefficient 
is very small at ranges near the effective range of both 
forces. For example, in Case I, with a force separation of 
1900 meters, the X force coefficient, or kill-rate, is given 
by 

3(1900) = Bo (1- r/R 0 ) 2 

= .005 (1- 1900/ 2000) 2 

= .0000125 . 
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The corresponding Y force kill-rate at this range is 

a (1900) = .001 . 

Thus, it is advantageous for the X force to close the range 
between forces as quickly as possible in order to reduce the 
time spent at ranges at which the Y force has such a distinct 
kill-rate advantage. 

Figure 10 shows that when m=2 and n=l, the X force was 
the superior force in both Case I and Case II, regardless 
of the attack velocity. However, as in the case with in = l 
or 2 and n=0, a reduced attack speed resulted in fewer X 
force casualties, while, at the same time, the Y force was 
annihilated. This is explained, as before, by the fact that 
during the early stages of the battle, the X force attrition- 
rate coefficient, together with the initial X force strength 
advantage, completely dominates the Y force coefficient and 
strength. Thus, it is more advantageous to attack at slower 
speeds . 
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IV. CONCLUSIONS 



From the previous discussion, it is clear that the 
effects of attack velocity on the outcome of Lanchester- 
type engagements depend on the parameter values of the 
particular engagement, and especially on the form of the 
attrition- rate coefficients of each force. When the at- 
trition-rate coefficient of the attacking force dominates 
that of the defender in the early stages of the battle, 
the attack should be conducted at a slow speed to take ad- 
vantage of the kill- rate superiority enjoyed by the attack- 
ing force. Conversely, when the defender dominates the 
action in the early stages, it is advantageous to attack at 
a faster speed. Thus, when the kill-rate of a force does 
not dominate the kill-rate of the opposing force at all 
ranges, it is most advantageous for that force to hold the 
opposing force at ranges where it is subjected to the maxi- 
mum kill-rate advantage. This implies the obvious result 
that the best tactic to employ is to keep the enemy in the 
"maximum killing zone." 

The cases in which the attack speed had the least effect, 
if any, on the outcome of the battle, were the ones in which 
one, or both, of the forces had a constant attrition-rate 
coefficient. It is intuitively appealing that there would 
be few practical cases in which either of the forces in- 
volved would have a constant kill-rate. Therefore, in the 
more practical cases discussed, it appears that the speed of 
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attack does influence the outcome of a battle. However, no 
axiomatic statement describing this effect can be given, 
since the effect depends on the characteristics and capa- 
bilities of the opposing forces. 

In an actual combat situation, the opposing forces 
might wel] have different effective ranges. An area for 
further investigation would be a study of the effects of 
attack velocity when forces do not have the same effective 
range. As previously stated, Taylor [6] has developed a 
solution to Bonder's model for the case when forces have 
different effective ranges and linear dependence of thp 
attrition-rate coefficients. Taylor states [6] that his 
solution can be extended to situations where other range 
dependencies of the coefficients exist. Such an investiga- 
tion using Taylor's results would undoubtedly lend greater 
insight into the effects of weapon-system mobility and 
attrition-rate dependencies. 
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V . SUMMARY 



Analytic solutions to variable -coefficient Lanchester- 

/ 

type equations were used to examine the effects of attack 
velocity on the outcome of engagements between two homo- 
geneous forces. Constant and range dependent (linear and 
quadratic) attrition-rate coefficients were considered. 

The terminal numerical strength for each force was deter- 
mined for various values of attack velocity. These termi- 
nal strengths were determined for the cases when the opposing 
forces had the same type of kill-rates, e.g., both linear, 
and when they had different types of kill-rates. 

When the opposing forces had the same type of kill- 
rates, the Lanchester square law applied. An attacking force 
with initial numerical superiority, but with firepower in- 
feriority, was able to achieve terminal numerical superior- 
ity by attacking at a high speed. This result is a consequence 
of the Lanchester square law, and the fact that an increase 
in attack velocity is equivalent to a decrease in the dura- 
tion of the battle. 

When the opposing forces had different types of kill- 
rates, the effects of attack velocity were dependent on the 
forms of the kill-rates. If the kill-rate of the attacking 
force dominated that of the defender, the attacker could 
conserve his force by attacking at slower speeds. Converse- 
ly, if the kill-rate of the defender was dominate, it was 
advantageous for the attacker to attack at higher speeds. 
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TABLE I . VALUES OF INPUT PARAMETERS 





Case I 


Case II 


X Force Initial Strength 


30 


30 


Y Force Initial Strength 


20 


20 


ao 


0.02 


0.02 


Bo 


0 . 005 


0.01 


Effective Range, R g 


2000m 


2000m 


m ,n 


0,1,2 


0,1,2 
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ATTRITION -RATE COEFFICIENTS ATTRITION-RATE COEFFICIENTS 





Figure 1. Attr it ion- rate coefficients, a(r) and 3(r), as 
a function of range for various values of m and 

n . 
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TERMINAL FORCE STRENGTH TERMINAL FORCE STRENGTH 



30 4 



Case I 



Figure 2 



10 - 




ATTACK VELOCITY (METERS PER SECOND) 
Case II 




ATTACK VELOCITY (METERS PER SECOND) 



Effect of attack velocity on terminal force 
strength; m=0, n=0. 
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TERMINAL FORCE STRENGTH TERMINAL FORCE STRENGTH 



Case I 




ATTACK VELOCITY (METERS PER SECOND) 



Case II 




ATTACK VELOCITY (METERS PER SECOND) 



Figure 3. Effect of attack velocity on terminal force 
strength; m=l, n=l. 
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TERMINAL FORCE STRENGTH TERMINAL FORCE STRENGTH 



Case I 




ATTACK VELOCITY (METERS PER SECOND) 



Case II 




ATTACK VELOCITY (METERS PER SECOND) 



Figure 4. Effect of attack velocity on terminal force 
strength; m=2, n-2. 
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TERMINAL FORCE STRENGTH TERMINAL FORCE STRENGTH 



Case I 




ATTACK VELOCITY (METERS PER SECOND) 



Figure 5. Effect of attack velocity on terminal force 
strength; m=0, n=l. 
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TERMINAL FORCE STRENGTH TERMINAL FORCE STRENGTH 



Case I 




ATTACK VELOCITY (METERS PER SECOND) 
Case II 




ATTACK VELOCITY (METERS PER SECOND) 



Figure 6. Effect of attack velocity on terminal force 
strength; m=0, n=2. 
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TERMINAL FORCE STRENGTH TERMINAL FORCE STRENGTH 



Case I 




ATTACK VELOCITY (METERS PER SECOND) 




ATTACK VELOCITY (METERS TER SECOND) 



Figure 7. Effect of attack velocity on terminal force 
strength; m=l, n=0. 
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TERMINAL FORCE STRENGTH TERMINAL FORCE STRENGTH 



Case I 




ATTACK VELOCITY (METERS PER SECOND) 
Case II 




ATTACK VELOCITY (METERS PER SECOND) 



Figure 8. Effect of attack velocity on terminal force 
strength; m=2, n=0. 
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TERMINAL FORCE STRENGTH TERMINAL FORCE STRENGTH 



Case I 



Figure 9 




ATTACK VELOCITY (METERS PER SECOND) 
Case II 




ATTACK VELOCITY (METERS PER SECOND) 



Effect of attack velocity on terminal force 
strength; m=l, n=2. 
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TERMINAL FORCE STRENGTH TERMINAL FORCE STRENGTH 



Case I 




ATTACK VELOCITY (METERS PER SECOND) 
Case II 




ATTACK VELOCITY (METERS PER SECOND) 



Figure 10. Effect of attack velocity on terminal force 
strength; m=2, n=l. 
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PROGRAM TO DETERMINE EFFECT OF ATTACK VELOCITY 



COMMON I PT, IOUT, ISTEP 
CATA H4END/4HEND / 

MAIN PROGRAM TO COMPUTE SERIES. BOTH WEAPON SYSTEMS HAVE THE 
SAME EFFECTIVE RANGE. 

10 FORMAT ( 214 ) 

20 FORMAT (A4,3X, 2 F6. 1,2F8.4,4I3,2F7.1,F9.1) 

30 FORMAT ( 1H1 ,5X, 1 LANCHESTER ATTRITION THEORY. - BOTH WEA 
1 PON SYSTEMS HAVE THE SAME EFFECTIVE RANGE . ' // /2H ,5X,' 
2 CASE NUMBER' , I 6//1H ,5X,7H X0=,F12.2/1H ,5X,7H Y 

3C= , F 12 . 2/1H ,5X,7HALPHA0=, F12.4/1H ,5X,7H BETA0=,F12.4 
4/1H , 5X , 7H MMIN=,7X, I5,5X,5HMMAX=,7X, I5/1H , 5X,7H NM 
5IN=,7X, 15. 5X,5HNMAX=,7X, I5/1H , 5X,7HDELTAR=, F12.2/1H , 
65X , 7 HR A NGE 0= » F 12.2/1H ,5X,7H V=,F12.2! 

40 FORM AT ( 1H1 ,5X, ' LANCHESTER ATTRITION PR0FILEV/1H ,9X,' 

1 BOTH WEAPON SYSTEMS HAVE THE SAME EFFECTIVE RANGE. '//I 

2 H ,8X, ' ALPHAO- * .F10.4/1H » 9X, ' BETAO= ' .F10.4/1H ,13X,'M 
3=' , 15/1 H ,13X, ' N=' , I5/1H , 13X , • V= ' , F8 . 2/ /// 1H ,5X,'IST 
4 E P ‘ ,4X,' RANGE' ,7X,'X(R)',6X,'Y(R)',7X,'F(X,Y)',10X,'IS 
5TEP 1 / > 

50 FORMAT ( 1H , 5 X , I 3 , 3 X , F 8 . 1 , F 1 1 .4 , F 1 0 . 4 , E 14 . 6 , 6 X , I 3 )' 

READ! 5 , 10 ) IPT, IOUT 
ICASE=0 

300 READ! I PT .20 ) FL AG, XO » YO , ALPHA3 » BET AO ,MM I N. MMAX , NMI N , NMA 
1X,DELTAR,V,RANGE0 
IF( F LAG.EQ.H4END ) GO TO 100 
I C A SE = I C AS E+ 1 

WRITE I IOUT ,30) I C AS E , XO , YO , ALPH AO , BET AO , MM I N , MMAX , NM I N , 
INMAX, DEL TAR, RANGEO, V 
DO 200 M=MM I N, MMAX 
DO 200 N=N M I N, NM AX 
R = RANG EO-DEL TAR 
I STFP= 1 

WRITE! I OUT, 40) A LP HAO , B E TAO , M , N , V 
400 CALL XYSER(XO,YO,ALPHAO,BETAO,M,N,R, RANGEO, V,X,Y) 
FXY~BETAO*X*X- ALPH AO*Y*Y 
WRITE! IOUT, 50) I ST E P , R , X , Y , F XY , I ST EP 

IF! .NOT. ( (X.GT.O.O) . AND . ( Y .GT .0 . 0 ) . AND . ( R . GT . 0 . 0 ) ) ) GO 
1T0 200 
R=R-DELTAR 
ISTFP=I STEP+1 
GO TO 400 
200 CONTINUE 
GO TO 300 
100 CONTINUE 
END 



SUBROUTINE XYS ER ( XO , YO , ALP HAO , BE TAO , M , N , R , RANGEO , V , X , Y 

1 ) 

THIS SUBROUTINE CALCULATES X ( R ) AND Y ( R ) ACCORDING TO THEIR 
SERIES SOLUTION. 

IMPLICIT REAL*8(C) 

COMMON IPT, IOUT, ISTEP 

10 FORMAT ! 1H ,5X, 'SERIES FOR X ( R > HAS NOT CONVERGED. NEXT 
1 ENTRY IS IN ERROR. ' ) 

20 FORMAT ( 1 H ,5X, 'SERIES FOR Y ( R ) HAS NCT CONVERGED. NEXT 
1 ENTRY IS IN ERROR. ' ) 

DATA EPS/1. 0E-06/ 

500 P=( M+l .0 )/ (M+N+2 .0 ) 

Q=l-P 

S = ( M+ N+2 . ) /2.0 

COMP 1= ( { RANGEO *SQRT ( ALP HAO* BE TAO ) ) / ( V*2 . *S ) ) **2 
C0MP2= 1 . 0-( R/RANGEO ) 

CONST 1=CCMP1*( COMP 2** (2*S ) ) 

C0NST2=C0MP1*( C0MP2** (2*P*S ) ) 

C0NST3 = CCMP1*( COMP 2** (2*Q*S) ) 

SOLUTION FOR X ( R ) 
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6C0 TERM 1=1. O/GAMMM 1.0-P) 

TERM2=(C0MP2**(2*P*S) ) / GAMMA ( 1 .0 + P) 

SUM1=TERM1 
SUM2=T ERM2 
DO 100 K=1 t 100 

TERM1=TERM1*C0NST1*GAMMA(K-P )/ ( K*GAMM A ( K+ 1 . 0- P ) ) 

TERM2 = TERM2*C0NST2*GAMMA(K + P )*<C0MP2**( 2*K*S) ) / ( K^GAMM 
1 A ( K+l .0+P)*C0MP2** (2*(K-1 ) *S+2*P*S) ) 

SUM 1 =SUM 1+TERM 1 
S UM2 = S U M 2 + T E RM 2 

I F < AMAX1( ABStTERMl ) ,ABS(TERM2) ) .LT.EPS) GO TO 200 
ICO CONTINUE 

WRIT E ( I OUT ,10) 

200 X=XG*GAMMA( 1 « 0 -P ) * SUM1 -R ANGEO* ALPHA0*Y0*GAMM A ( 1 .0+P )*S 
1 U M 2 / ( V* ( M+l .0) ) 

SOLUTION FOR Y ( R ) 

TERM3= 1 . O/GAMM A ( 1 ,0-Q ) 

TERM4=CCMP2**(2*Q*S) /GAMMA ( 1 .0 + Q) 

SUM3=TERM3 
SUM4=T E RM4 
DC 300 J=1 , 100 

TERM3=TERM3*C0NST1*GAMMA( J-Q)/ ( J*GAMMA( J + l .0-Q) ) 
TERM4=TERM4*C0NST3*GAMMA( J+Q)*C0MP2** (2*J*S) / ( J*GAMMA( 
1 J + l. 0+Q »*CCMP2** (2*( J-l ) *S+2*Q*S) ) 

SUM3=SUM3 + TERM 3 
SUM4=SUM4+TERM4 

I F ( AMAX1 ( A3S (TERM3 ) , ABS< TERM4) ) ,LT .EPS ) GO TO 400 
3C0 CONTINUE 

WRITE! I OUT ,20) 

400 Y=YO*GAMMA< 1 . 0-Q ) * SUM3-RANGE0* B ET AO*XO*G AMM A ( 1 .0+Q )*SU 
1M4/(V*(N+1.0)) 

RETURN 

END 
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